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Using simple dimensional arguments for both spiral and elliptical galaxies, we 
present formulas to derive an estimate of the halo spin parameter A for any real 
galaxy, in terms of common observational parameters. This allows a rough estimate 
of A, which we apply to a large volume limited sample of galaxies taken from the 
SDSS data base. The large numbers involved (11,597) allow the derivation of reliable 
A distributions, as signal adds up significantly in spite of the errors in the inferences 
for particular galaxies. We find that if the observed distribution of A is modeled with 
a log-normal function, as often done for this distribution in dark matter halos that 
appear in cosmological simulations, we obtain parameters Ao = 0.04 ± 0.005 and 
o\ = 0.51 ± 0.05, interestingly consistent with values derived from simulations. For 
spirals, we find a good correlation between empirical values of A and visually assigned 
Hubble types, highlighting the potential of this physical parameter as an objective 
classification tool. 

Key words: galaxies: statistics - galaxies: formation - galaxies: fundamental param- 
eters - galaxies: structure - cosmology: observations - cosmology: miscellaneous 



1 INTRODUCTION 

Perhaps one of the most useful tools for describing the na- 
ture of observed galactic populations is the luminosity func- 
tion. The relative numbers of galaxies of different luminosi- 
ties offer a quantitative description of the result of the struc- 
ture formation scenario of the universe, nowadays available 
for inspection even as a function of redshift. The luminos- 
ity function, or mass function of galaxies, once a modeling of 
star formation histories and dust is included to yield mass to 
light ratios, is indeed one of the principal constraints against 
which cosmological simulations and structure formation the- 
ories in general are tested. Often star formation efficiencies 
and histories are calibrated through a fitting of predicted 
halo mass functions and observed luminosity functions e.g. 
Davis et al. (1985) and references thereof. 

Galaxies however, have many more properties than just 
their mass, most notably, galactic classification schemes cen- 
ter on sorting galaxies as ellipticals, lenticulars and spirals, 
subdividing the latter again into types, typically along the 
Hubble sequence. In spite of the many well known virtues of 
this classification scheme, its somewhat subjective, qualita- 
tive and relative nature, has made it difficult to use in com- 
parisons against cosmological simulations, where the 'type' 



of a modeled galaxy is rather difficult to assess, in terms of 
Hubble's classification scheme. 

Inspired by varied theoretical studies which invariably 
identify the A spin parameter of a host halo as the principal 
physical parameter in determining the morphological and vi- 
sual characteristics of a spiral galaxy, which are then subjec- 
tively integrated into the qualitative assignment of 'type', in 
Hernandez & Cervantes-Sodi (2006) two of us derived a sim- 
ple estimate of A for any observed spiral. There it was shown 
that the scalings of the derived A against various type deter- 
mining properties such as colour, disk thickness and bulge to 
disk ratio, are comparable to the corresponding scalings of 
these parameters against Hubble type, using a large sample 
of nearby spirals. 

It is interesting that a generic prediction of cosmolog- 
ical N-body simulations over the past 2 decades has been 
the functional form and parameters which define the pre- 
dicted distribution of A for dark matter halos e.g. Bullock 
et al. (2001). Nevertheless, it has not been easy to test this 
prediction directly, as A is not a straight forward observ- 
able feature of a real galaxy. Not only is a measurable A 
distribution relevant as a test of the general structure for- 
mation scheme, but also as a further way of independently 



2 X. Hernandez, Changbom Park, B. Cervantes- Sodi and Yun- Young Choi 



constraining cosmological parameters, which to a certain ex- 
tent alter its details e.g. Gardner (2001). 

Having the approximate estimates of Hernandez & 
Cervantes-Sodi (2006) (henceforth HC06), here we set out 
to measure empirical distributions of galactic A spin param- 
eters. Firstly, given the approximate nature of these esti- 
mates, the only way of having signal adding up to some- 
thing significant is through the use of very large samples. 
Also, if we want a meaningful comparison against cosmolog- 
ical models, use of a volume limited sample is crucial. These 
two constraints drive us inevitably to the Sloan Digital Sky 
Survey (SDSS). Two of us have worked extensively with this 
database, and recently constructed an interesting morphol- 
ogy segregation scheme of galaxies in the SDSS in Park & 
Choi (2005), henceforth PC05. Here we use a volume limited 
sample from the SDSS having galaxies with redshifts in the 
interval 0.025 < z < 0.055, corresponding to distances of be- 
tween 75fr~ 1 Mpc and 162.57/iT 1 Mpc, assuming a WMAP 
cosmology of Qm = 0.27, £Ia = 0.73 and h — 0.71, which we 
keep throughout. 

The following section includes a review of the derivation 
of A spin parameters for spirals of HS06, and gives our results 
for the spirals in our SDSS sample. The correspondence be- 
tween A for spirals, galactic type through visual inspection 
and the colour and colour gradient morphology segregation 
scheme of PC05 is also made in Section 2, where we also 
present empirical distributions of galactic A spin parameters 
for spirals. These are well fitted by a log-normal function, 
and give parameters in accordance with recent cosmological 
simulations. In section 3 we extend the dimensional estimate 
of A for ellipticals, and construct and analyze the empirical 
distributions of A for the complete sample. Our results are 
summarized in 4. 



2 A DISTRIBUTIONS FOR SPIRAL GALAXIES 

We are concerned with the determination of the dimension- 
less angular momentum parameter of a galactic halo: 

where E, M and L are the total energy, mass and an- 
gular momentum of the configuration, respectively. None of 
the quantities which appear in the above definition are ac- 
cessible to direct observation. In HC05 we derived a simple 
estimate of A for disk galaxies in terms of observational pa- 
rameters, which we summarize below. 

We shall model only two galactic components, the first 
one a disk having a surface mass density E(r) satisfying: 



E(r) = E e- r/ ^, 



(2) 



Where r is a radial coordinate and So and R d are two 
constants which are allowed to vary from galaxy to galaxy. 
The total disk mass is now: 



M d = 2-KT. R 2 d . 



(3) 



The second component will be a fixed dark matter halo 
having an isothermal p(r) density profile, and responsible for 
establishing a rigorously flat rotation curve Vd throughout 
the entire galaxy, an approximation sometimes used in sim- 
ple galactic evolution models e.g. Naab & Ostriker (2005). 



Under a flat rotation curve the disk specific angular momen- 
tum from eq(2) will be Id = 2VdRd, and the halo density will 
satisfy: 



p(r) 
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(4) 



having a halo mass profile M(r) oc r. We define the disk 
mass fraction as F — Md/Mn, expected to be of order 1/10 
or smaller (e.g. Flores et al. 1993, Hernandez & Gilmore 
1998), we shall use global parameters for the entire system 
indistinctly from halo parameters, consistent with having 
ignored disk self gravity in eq(4). 

We must now express A in terms of structural galac- 
tic parameters readily accessible to observations. First we 
assume that the total potential energy of the galaxy is dom- 
inated by that of the halo, and that this is a virialized grav- 
itational structure, which allows to replace E in equation(l) 
for W/2, one half the gravitational potential energy of the 
halo, given by: 



W ■ 



-V d 'M H . 



(5) 



Assuming that the specific angular momenta of disk and 
halo are equal, Id = Ih (as would be the case for an ini- 
tially well mixed protogalactic state, and generally assumed 
in galactic formation models e.g. Fall & Efstathiou 1980, Mo 
et al. 1998) we can replace L for Mnld- Introducing this last 
result together with eq(5) into eq(l) yields: 



GM H ' 



(6) 



Finally, we can replace Mh for M d /F, and introduce a 
disk Tully-Fisher relation: 



M d = A TF vr 

into eq(6) to yield: 



A = 



2 1 ' 2 F 
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(7) 



(8) 



The existence of a general baryonic Tully-Fisher rela- 
tion between the total baryonic component and Vd of the 
type used here, rather than a specific Tully-Fisher relation 
involving total magnitude in any particular band, is sup- 
ported by recent studies, in general in agreement with the 
3.5 exponent we assume (e.g. Gurovich et al. 2004 or Kregel 
et al. 2005 who find 3.33 ± 0.37). All that remains is to 
choose numeric values for F and Atf, to obtain an esti- 
mate of the spin parameter of an observed galaxy in terms 
of structural parameters readily accessible to observations, 
R d and V d - 

Taking the Milky Way as a representative example, 
we can use a total baryonic mass of 1 x 10 n Mo and 
M H = 2.5 x 1O 12 M (where the estimate of the bary- 
onic mass includes all such components, disk, bulge, stel- 
lar halo etc., e.g. Wilkinson & Evans 1999, Hernandez et 
al. 2001) as suitable estimates to obtain F = 1/25. For a 
rotation velocity of Vd = 220km/ s, the above numbers im- 
ply Atf = 633Mo(fcm/.s) -3 ' 5 , this in turn allows to express 
eq(8) as: 



A = 21.8 



Rd/kpc 



(Vd/kms- 1 )*/ 2 ' 



(9) 
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For the Galactic values used above, equation (9) yields 
Amw = 0.0234, in excellent agreement with recent esti- 
mates of this parameter, e.g. through a detailed modeling 
of the Milky Way within a CDM framework Hernandez et 
al. (2001) find X M w = 0.02. 

First, we note that by construction, what we are esti- 
mating is not strictly A, but what has been defined as A', 
the equivalent A for a singular truncated isothermal halo. 
The relation between A and A' is slightly a function of halo 
structure, for example, for NFW profiles with concentration 
of 10, one gets A' = 0.9A e.g. Mo ct al. (1998). This slight 
difference is smaller than the error introduced through the 
dispersion in the Tully-Fisher relation used, and will not 
be considered further, in any case, it is more rigorously A' 
and not A that we will be talking about. Also, it is clear 
that eq(9) is at best a first order approximation, and not a 
precise evaluation of A for a real galaxy. 

However, in HC06 it was shown that this A parameter 
shows a one-to-one correlation, with small dispersion, when 
compared to the input A of detailed galactic formation mod- 
els from two distinct cosmological groups. Also, the scalings 
seen against colour, disk thickness and bulge to disk ratios, 
for a sample of nearby spirals from Courteau (1996, 1997), 
de Grijs (1998) and Kregel et al. (2002), similar to what is 
seen against Hubble type, highlights the use of this parame- 
ter as a physical classification scheme. This is not surprising, 
as it is A as defined by eq.(l) what has been repeatedly iden- 
tified in analytic and numerical studies of galactic formation 
as the principal determinant of galactic type e.g. Fall & Ef- 
stathiou (1980), Flores et al. (1993), Firmani et al. (1996), 
Dalcanton et al. (1997), Zhang & Wyse (2000), Silk (2001), 
Kregel et al. (2005). 

For example, the ratio of disk scale height to disk scale 
length, h/Rd, is one of the type defining characteristics of a 
galaxy which it is easy to show, will scale with A. Starting 
from Toomre's stability criterion: 



Q{r) = 



n(r)a(r) 
7rGS(r) ' 



(10) 



were /t(r) and a(r) are the epicycle frequency and veloc- 
ity dispersion at a given point in the disk, assuming a thin 
disk, virial equilibrium in the vertical direction (Binney and 
Tremaine 1994) yields a relation between h and E, 



h = 



2tyGT,' 



(11) 



We use this relation for h to replace the gas velocity 
dispersion appearing in equation (10) for a combination of 
h and the surface density. The dependence on E is replaced 
by one on Md and Rd through the disk profile. Replacing 
n(r) for \/2Vd/r, and Vd for A, Rd and Mr through eq(6) 
to get a new expression for the Toomre's stability criterion, 
which evaluating radial dependences at r = R d yields, 



(12) 



With F = 1/25, evaluating at Q = 1, the stability 
threshold suggested by self-regulated star formation cycles, 
(e.g. Dopita & Ryder 1994, Koeppen et al. 1995 and Silk 
2001) the ratio h/Rd is obtained as: 



1 



R d 390A ' 



(13) 



a simplified version of the result of van der Kruit (1987) . 
For the Galactic value derived above of Xmw = 0.0234, 
equation (13) gives Rd = 9h, not in conflict with parame- 
ters for the Milky Way. For galaxies with large values of A 
we expect thin systems, while galaxies with small values of 
A will show thick disks, as the observed trend of (h/Rd) de- 
creasing in going from early-type disks to late type galaxies 
e.g. de Grijs (1998), Yoachim & Dalcanton (2006). 

Having a way of estimating A, we now need to define 
the optimal sample to use. As already mentioned, requir- 
ing as large as possible a volume limited sample makes the 
SDSS an ideal database. Since most published cosmological 
distributions of A report data at z = 0, we will start by 
taking a volume-limited, absolute magnitude limited sam- 
ple with M r — 5 log h < —18.5, in a low redshift range, 
0.025 < z < 0.055, allowing an optimal comparison with 
reported cosmological studies. This sample contains 31,685 
galaxies for which exponential disk scales, absolute magni- 
tudes, velocity dispersions, de Vacouleurs radii and eccen- 
tricities have been determined (Choi et al. 2006). It is im- 
portant to note that the above structural parameters were 
determined in red and IR bands, in the spirit of obtaining 
global mass distribution parameters, it is clear that this is 
what is required for the measure of Rd necessary in eq(9). 

The first step is to identify the spiral galaxies in our 
sample. This is done using the colour, colour-gradient and 
concentration criteria developed in PC05, where extensive 
testing and corroboration of the morphological segregation 
criteria were preformed against large training sets of visually 
classified galaxies. This procedure yields 21,184 spirals and 
10,501 ellipticals. 

Since we have no rotation curves for the spiral galax- 
ies, as required by equation(9), we must infer this through 
the absolute magnitudes of the observed systems, and use 
of an appropriate Tully-Fisher relation. Hence, internal ab- 
sorption in edge-on disks is a problem we must avoid. We 
prune the original sample to leave only spiral galaxies hav- 
ing axis ratios > 0.6, this ensures only relatively face-on 
disks remain, minimizing errors in the absolute magnitudes 
used (see Choi et al. 2006 for the choice of axis ratio cut). 
Also, we use the red band Tully-Fisher relations of Barton et 
al. (2001) to assign rotation velocities to observed galaxies. 
This relation shows a good fit to the data only in the range 
—20 > Mr > —22.5, so we must limit our galaxies to those 
falling in this range. The corresponding ranges we impose on 
inferred rotation velocities are 80 < Vr < 430 in km/s, well 
within the range of applicability of the Tully-Fisher relation 
we are using. Still, the dispersion in this Tully-Fisher, plus 
that in the baryonic Tully-Fisher used in deriving eq.(9), 
leave us with a 25% uncertainty in our individual estimates 
of A for spirals. 

After applying the two cuts described above, we are left 
with 0.366 of the original 21,184 spirals, 7,753 disk systems, 
still large enough to yield significant statistical information 
on the distribution, and small enough to include only galax- 
ies for which our estimates of A are most accurate. Use of 
eq.(9) for this sample yields a collection of 7,753 values of A 
for a complete volume limited SDSS sample. 

Having now a set of inferred, observational values of A, 
we first explore the correlations of A for spirals only, against 
the colour and colour-gradient information found by PC05 
to be an accurate indicator of galactic type. Figure (1), 
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Figure 1. Some of the galaxies in our sample on a colour, colour- 
gradient plane, with photos of galaxies representative of those 
found in each region, as visually classified by Park & Choi (2005). 
Red dots are early type galaxies and others are late types. 



taken from PC05, shows some of the galaxies in our sam- 
ple, with photos of galaxies representative of those found in 
each region, as visually classified by PC05. Notice the well 
denned cluster of ellipticals centered on A(g — i) = —0.04, 
u — r = 2.82. Figure (2) shows all our spirals in a colour, 
colour-gradient plane, with the shading giving the average 
values of A within each shaded square. We see that disks 
with high values of A are found in the lower left hand area, 
whilst disks having low values of A populate the right and 
upper regions. It is interesting that PC05 find precisely the 
same segregation pattern for late and early spirals, respec- 
tively, reinforcing the results of HC06 that the quantitative 
and objective A parameter constructed in eq.(9) is a good 
physical classification parameter for spirals, reproducing the 
broad trends of the classical subjective, qualitative Hubble 
sequence. 

Next we use our collection of values of A from our 
pruned sample of face-on spirals to construct a histogram, 
shown in figure (3) by the broken curve. The shape of this 
histogram suggests a log-normal distribution: 



P(A ,a A ;A)dA = 



exp 



ln 2 (\/\ Q ) 



2*1 



dX 
T 



(14) 



The moments of the above distribution are analytical, 
which allows one to express: 

f n °° P(A)AdA 



J °°P(A)dX 
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Jo P ( A ) dA 



from which one solves: 
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Figure 2. Spiral galaxies in our sample on a colour colour- 
gradient plane, the shading shows the average values of A in each 
shaded region, using eq. (9). Regions of high and low values of A 
closely correspond to regions populated by late type and early 
type spirals, respectively, as visually classified by Park & Choi 
(2005). 



Equating the 1st and 2nd order moments of our em- 
pirical distribution to Mi and M2, we find the best fit log- 
normal curve to our distribution for spiral galaxies. This is 
given by the smooth curve in figure (3), where the empirical 
histogram has been rescaled to a unit integral. The parame- 
ters of the best fit distribution are: Ao = 0.0585, a\ = 0.446. 
Both the form of this distribution and the particular values 
for the parameters we find are in good agreement with re- 
sults from cosmological simulations, e.g. Shaw et. al (2006). 

From figure (3) it is apparent that the log-normal func- 
tional form is a good description of our inferred data, a 
more formal confirmation is provided by a K-S test. We 
construct the cumulative distribution functions for both the 
log-normal P(A) and our data out to A = 0.2, shown in fig- 
ure (4) by the solid and dotted curves, respectively. These 
yield a maximum difference D = 0.029, and a significance 
given by: 



Qks(X) =2j2(~ 1 Y' 1 ^P(^j 2 X 2 



(18) 



j=i 



where X — N X ^ 2 D, N the total number of points in 
our sample, 7,746, giving X = 2.55 and a significance level 
of Qks = 0.000004. We hence see that our inferred distri- 
bution of galactic A's is consistent with having come from 
the best fit log-normal -P(A), to a very high degree. We note 
the important precedent of Syer et. al (1999) who using sim- 
ilar estimates for A for spiral galaxies, using a sample from 
Lauberts (1982), find a distribution of values of A well fitted 
by a log-normal function, with parameters Ao = 0.05 and 
a\ = 0.36. The small differences with our inference might 
come from their lack of a large volume limited unbiased sam- 
ple, their lack of an attempt at modeling A for ellipticals, 



Empirical distributions of X 5 




0.05 0.1 0.15 0.2 0.25 

A 



Figure 3. Distribution of values of A for 7,753 spirals using eq.(9), 
binned into 150 intervals, broken curve. Best log-normal fit to the 
data, having parameters Ao = 0.0585 and u\ = 0.446, solid curve. 

which further limited the scope of comparisons which could 
then be made against cosmological n-body models, or their 
use of a disk formation model calibrated directly from cos- 
mological n-body simulations. 

From equation(8) we see that our inferred A's are pro- 
portional to the constants (F/Atf), since the validity and 
constancy of the T-F relation is well established, we can 
think of our inferred A's as A(25P), where F is the baryon 
fraction of a spiral system. We have shown that under the 
simplest assumption of a constant F = 1/25, a log-normal 
distribution results from the data, clearly, any other con- 
stant F will also yield a log-normal distribution, with the 
same cr\, and a mean of 25FAo- Also, notice that most cos- 
mological simulations yield distributions of A being indepen- 
dent of mass (e.g. Bullock et al. 2001) which if true, ensures 
our inferences will not be skewed by the velocity cuts applied 
to the spirals. 

A baryon fraction being a strong function of Rd or Vd 
however, could easily result in a best fit P(A) no longer 
looking like a log-normal distribution. In the absence of any 
definitive handle on F or its possible variations, we limit 
ourselves to concluding that a constant F implies a log- 
normal P(A) for the spirals from the SDSS we studied, while 
varying F{Rd,Vd) could doubtlessly imply something else. 



3 A DISTRIBUTIONS FOR THE COMPLETE 
SAMPLE 

Before attempting an empirical total distribution of A's we 
must derive an estimate for this parameter for an ellipti- 
cal galaxy, we propose a model equivalent to the one used 
for spirals, including only two components; a baryonic mat- 
ter distribution with a Hernquist density profile and a dark 
matter halo having a singular, truncated isothermal profile. 
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Figure 4. Cumulative distribution functions for both our spiral 
galaxies data and the best fit log-normal P(A), dotted and solid 
curves, respectively. 



In order to calculate the energy of the galaxy, we as- 
sume that the total potential is dominated by the energy 
of the virialized dark matter halo. In principle, we expect 
that halos were elliptical galaxies are found are no different 
from those of disk galaxies (White & Rees, 1978, Kashlin- 
sky, 1982). In the case of a disk galaxy, the halo density 
profile has the form AnGp(r) = (V c /r) 2 , with potential en- 
ergy —V^Mh- The velocity V c is the circular velocity of an 
element in centrifugal equilibrium with the halo. In an ellip- 
tical, nothing is actually moving at that velocity, which still 
can be used to characterize the halo. The dependence on ve- 
locity can be changed to one on mass using again a baryonic 
Tully-Fisher relation (Gurovich et al, 2004): M b = AtfV^ 5 , 
thinking that the dynamics of a halo host of an elliptical 
galaxy are the same as those of disk galaxies. Finally, we 
define the baryonic mass fraction as F — M^/Mh to ex- 
press the mass of the halo in terms of the baryonic mass. 
Validating the above assumptions, we note the interesting 
new results of Treu et. al (2006) and Koopmans et. al (2006) 
who find trough dynamical studies of velocity dispersion in 
ellipticals, and direct inferences on the halos of the studied 
ellipticals through lensing measurements, 1) the presence of 
extensive dark halos and 2) that these halos have equiva- 
lent velocity dispersions equal to those of their galaxies, the 
dynamical equivalent of a flat rotation curve for ellipticals. 

For the angular momentum, if we assume that the spe- 
cific angular momenta of dark matter and baryons are equal, 
h = Ih, we can obtain the angular momentum of the entire 
configuration using only the directly observable component. 
It is important to remember that the baryonic component 
is susceptible of dissipating, while dark matter is not, this 
will affect the above assumption and our A estimates, be- 
ing the value obtained a lower limit. In order to obtain the 
angular momentum of the baryonic distribution, we need to 
know the rotation velocity, information not always available 



6 X. Hernandez, Changbom Park, B. Cervantes- Sodi and Yun-Young Choi 



due to technical observational difficulties and the fact that 
in many elliptical galaxies, the rotational velocity is at best 
of the same order as the velocity dispersion (Pinkney et al. 
2003, Sarzi et al. 2006). These forces us to determine the an- 
gular momentum of the system using some other observable 
parameter. In many works (e.g. Binney 1978, Franx et al. 
1991) it is shown that the ellipticity of a galaxy is related to 
its rotation velocity, more precisely to its angular momen- 
tum. By dimensional analysis we expect the specific angular 
momentum of the galaxy be proportional to (GM^a) 1 / 2 and 
to the eccentricity e of the galaxy in the i-band, as shown by 
Gott & Thuan (1976) and in Marchant & Shapiro (1977). In- 
troducing a proportionality constant C related to the details 
of the matter distribution, the specific angular momentum 
is given by: 



l b = CeVGM b a, 



(19) 



with a the de Vacouleurs radius of the galaxy. The constant 
C can be calibrated from cases where we know the actual 
angular momentum of the system. For a small sample of 22 
closely studied elliptical galaxies, from Halliday et al. (2001) 
and Pinkney et al. (2003) we calculated the actual angular 
momentum, for a Hernquist profile and using the reported 
values of rotational velocity. Comparing with the angular 
momentum given by eq.(19), we obtain C = 0.114. 

For the galactic mass, using again the assumption that 
the gravitational structure is virialized, we compute the dy- 
namical mass at R$q, the half light radius, using the expres- 
sion 



dyn 



(1.65a) 2 R 50 
G 



(20) 



from Padmanabhan et al. (2004) , where a is the velocity 
dispersion. For elliptical galaxies it is well known that the 
presence of dark matter in the inner part of the galaxy is 
negligible so the total baryonic mass will be simply 2Md yn - 

Now we have the tools to estimate the energy, the an- 
gular momentum and the mass of an elliptical galaxy. If we 
introduce this information into equation 1, we obtain: 



A = 



0.051e(a/fcpc) 2/7 
(a/kms- 1 ) 3 / 7 ' 



(21) 



Where we used values for Atf and F as calibrated in 
HC06. The above is the equivalent expression of eq.(9) for 
elliptical galaxies, and has errors of 30%, due to the disper- 
sion in the observational relations used. 

We must take into account that the values for eccentric- 
ity e and radius a, required must be the intrinsic parameters 
of the system, not the observed ones, as these will be system- 
atically smaller due to projection effects. Also in many ellip- 
ticals, shape is the result not exclusively of rotational sup- 
port, but also of dynamical pressure support (i.e. anisotropic 
velocity dispersion), with varying degrees of relative impor- 
tance. Hence, unknown dissipational and projection effects 
make eq.(21) a lower limit, whilst the participation of pres- 
sure support make it an upper one. To what extent the above 
compensate each other is unclear, we shall therefore assume 
eq.(21) must be multiplied by a factor which remains a free 
parameter of our method, Xe- 

Since ellipticals will generally have much lower values of 
A than spirals, the relative mix of spirals and ellipticals will 
be an important factor in determining our final estimated 
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Figure 5. Specific angular momentum for elliptical galaxies 
from Pinkney et. al (2003) (triangles) and Halliday et. al (2001) 
(rhombi), determined directly from their measurements of rota- 
tion velocity and deprojected light profiles, x axis. The y axis 
shows values for the same quantity, but estimated through equa- 
tion (19) using C = 1, the solid line gives the best fit linear rela- 
tion. We then estimate C by forcing the best fit linear relation to 
coincide with the identity. 



distribution of galactic A's. We must hence leave this rela- 
tive fraction equal to that of the original sample, before the 
velocity and projection cuts were applied to the spirals, we 
further remove randomly selected ellipticals to reduce the 
total number of galaxies to 0.366 of the original number, 
leaving only 3,843 ellipticals. Our final sample is composed 
of a total of 11,597 galaxies, large enough to yield significant 
statistical information on the distribution. 

Using eq.(21) we estimate A for the ellipticals in the 
sample, after which we can add the two distributions to ob- 
tain the total empirical -P(A). As discussed previously, for 
ellipticals, a degree of uncertainty remains, in connection to 
the (varying) degree of pressure support, the degree of dis- 
sipation of angular momentum in the baryonic component, 
and projection effects, which to Oth order we have incorpo- 
rated as an unknown multiplicative factor Xe- For Xe = 1 
the result of adding the distribution of A parameters for spi- 
rals and ellipticals is a distinctly double peaked -P(A), with 
the ellipticals coming in as a gaussian distribution centered 
at A = 0.01. In the interest of comparing with simulations, 
we shall take as a working hypothesis that the full -P(A) 
should also be a log-normal distribution. We then adjust 
Xe so as to obtain a log-normal total distribution, with the 
restriction that the values of A for ellipticals should remain 
lower than those of spirals. The resulting histogram, requir- 
ing Xe = 2.3 is shown by the broken curve in figure(6), 
where the smooth solid curve gives a best fit log-normal 
curve to the full distribution, having parameters Ao = 0.046 
and cr\ = 0.535. 

Given that we have used empirical scaling relations, 
which show a measure of dispersion, we now take the pes- 
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simistic view that such dispersions are intrinsic to the data 
(pessimistic in the sense that a measure of this dispersion 
is undoubtedly observational noise), and hence, propagate 
as errors into our estimates of A. We now estimate what is 
the best underlying log-normal distribution, assuming that 
the histograms we obtained are only a degraded version of 
an intrinsic log-normal distribution. We determine the best 
intrinsic log-normal function through a full maximum likeli- 
hood analysis of the data. This is done for purposes of com- 
parison against predicted A distributions, we fix the func- 
tional form of the distribution, and set out only to evaluate 
the best fit parameters Ao and ax from eq.(14). 

We construct the likelihood function as the probability 
that an inferred set of n empirical values of A, {Ai(A)} n , 
might arise from a given model, (Xo,ax), as: 

71 poo 

C{\ ,oy,{\ t (\)} n ) = ^ln j P(A ,a A ;A)xA 1 (A)dA(22) 
i=0 Jo 

where: 



Ai(A) = 



F 



ai\/2ir 



exp 



(A-ApQ 2 
2a? 



(23) 



Aoj is the center of a density distribution, the nom- 
inal inferred A for the ith galaxy, at the assigned error, 
and Fi a normalization factor taking into account that no 
matter what the error, the domain of Xi must always be 
the positive values of A. Therefore, we normalize Ai(A) as 
a truncated Gaussian having unit integral, giving Fi = 
2/[l + erf (Aio/Va*)]. 

In the limit of the error in our inferred A tending to zero, 
Fi tends to 1, and A^(A) tends to a Dirac delta function, re- 
ducing the integral in equation(22) to an evaluation of P(X) 
at the inferred value. The advantage of a full likelihood for- 
mulation is that parameter inference can be preformed in a 
way which naturally incorporates the errors in the data sam- 
ple, without the need of binning the data, a process which 
intrinsically reduces the information content of the sample. 

Equation(22) is then evaluated over a fine grid of values 
of (Ao, ax), and the point where the maximum is found 
selected as the best fit model. For the 7,753 spirals the result 
is (Ao = 0.0517, ax ~ 0.362), with very small confidence 
intervals of ±0.0003 and ±0.004, respectively. Although our 
errors in individual A estimates are of 25% of A, it is thanks 
to having a sample running into the thousands that we can 
retrieve details of the A distribution with accuracy. 

Finally, we repeat the evaluation of eq.(22), but using 
this time the complete sample, of 11,597 spirals plus ellipti- 
cals. Here A for the ellipticals is taken as 2.3 x the estimate 
of eq.(21). This factor can not be much larger than what we 
are using, since then ellipticals would start overlapping sig- 
nificantly with spirals in their A distributions, as mentioned 
previously, if this factor is much reduced, the total A distri- 
bution becomes double peaked, with ellipticals appearing as 
a distinct population at very low A, which would be hard to 
explain. 

The results this time are Ao = 0.0394 and ax = 0.509. 
A comparison of the maximum likelihood model and the 
collection of inferred As is given in fig. (6), which is analogous 
to fig. (3), but includes the complete sample, binned into 150 
discrete intervals. A direct fit of eq.(14) to the full data gives 
A = 0.046, ax = 0.535. 




Figure 6. Distribution of values of A for 11,597 spirals plus el- 
lipticals using eq.(9) and eq.(21), binned into 150 intervals, bro- 
ken curve. Best log-normal fit to the data, having parameters 
Ao = 0.046 and ax = 0.535, smooth solid curve. Maximum like- 
lihood underlying log-normal distribution, assuming a 30% error 
in our inferred A's, with parameters Ao = 0.0394 and ax = 0.509, 
dotted curve. 



The formal errors in the method are again of the or- 
der of what was found for the spirals, but this time we are 
dominated by the unknown correction factor between our 
estimate and the actual A for ellipticals. This uncertainty, 
although bounded, dominates our final error estimates and 
yields ±0.005 in Ao and ±0.05 in ax- Another possible source 
of error in our estimates would be the existence of a large 
population of low surface brightness galaxies of high A, which 
would lead to larger values of both Ao and ax ■ 

The above results qualitatively agree with generic pre- 
dictions of structure formation models regarding the pre- 
dicted values for Ao and ax, e.g. Shaw et al. (2006) review 
recent results from the literature giving values in the range 
0.03 < A < 0.05 and 0.48 < a x < 0.64. The availability 
of empirical measurements of these parameters could serve 
as a further independent guideline for models of structure 
formation. 

We stress that for the full sample, the log-normal na- 
ture of both the inferred and the underlying A distributions 
is an assumption which we can not confirm, taken only to 
permit a comparison with models. Under this assumption, 
the parameters of the distribution can be derived from the 
data. An inference of the actual functional form of the to- 
tal distribution would require more solid constraints on the 
hypothesis introduced in the derivation of the A's for ellipti- 
cals. However, we do note that any similar scheme modeling 
ellipticals as residing in halos equivalent to those of spirals, 
and having significantly lower angular momentum than spi- 
rals, will also yield a P(X) not very different from what is 
shown in figure (6). 
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4 CONCLUSIONS 

We apply simple estimates of the A spin parameter to a large 
volume limited sample from the SDSS. The sample is split 
into ellipticals and spirals using colour, concentration and 
colour gradients. 

We find that for spiral galaxies, the average value of the 
inferred A correlates well with standard type, as determined 
by visual inspection. This highlights the potential of A as 
an objective and quantitative classification tool for spiral 
galaxies. 

For the spiral galaxies sample we obtain a distribution 
which is statistically consistent with a log-normal distribu- 
tion, as what is commonly found in cosmological n-body 
simulations, we find parameters Ao = 0.0585, a\ = 0.446. 

Ellipticals have, as expected, average values of A an or- 
der of magnitude lower than spirals. The details of their A 
distribution are harder to quantify, but average values of 
around 0.005 are derived. 

If the distribution of A parameters for the full sample 
is fitted by log-normal distribution, of the type found to 
reproduce the corresponding distribution of modeled halos 
arising in cosmological simulations, the parameters we find 
are: A = 0.0394 ±0.005 and a x = 0.509 ±0.05, derived this 
time from a sample of real galaxies. 
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